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We use particle entanglement spectra to characterize bosonic quantum Hall states on lattices, 
motivated by recent studies of bosonic atoms on optical lattices. Unlike for the related problem of 
fractional Chern insulators, very good trial wavef unctions are known for fractional quantum Hall 
states on lattices. We focus on the entanglement spectra for the Laughlin state at v = 1/2 for the 
non-Abelian Moore- Read state at v = 1. We undertake a comparative study of these trial states 
to the corresponding groundstates of repulsive two-body or three-body contact interactions on the 
lattice. The magnitude of the entanglement gap is studied as a function of the interaction strength 
on the lattice, giving insights into the nature of Landau-level mixing. In addition, we compare the 
performance of the entanglement gap and overlaps with trial wavefunctions as possible indicators 
for the topological order in the system. We discuss how the entanglement spectra allow to detect 
competing phases such as a Bose-Einstein condensate. 
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I. INTRODUCTION 

Optical lattices present unique opportunities to simu- 
late the physics of charged particles in strong magnetic 
fields While early proposals for artificial gauge fields re- 
lied on rotation to mimic the action of Lorentz forces by 
the Coriolis forces acting in the rotating frame opti- 
cal lattices provide a robust experimental setting where 
the effect of fields can be simulated by imprinting com- 
plex phases on tunnelling elements between neighbouring 
sitespjED xhe most elegant schemes for generating syn- 
thetic gauge fields in atomic gases rely on the use of 
Berry-phases resulting from a set of internal states sub- 
ject to a spatially varying optical dressing!^^ In particu- 
lar, these schemes enable the generation of high densities 
of flux per plaquette on the underlying lattice, giving ac- 
cess to a regime of strong correlation where exotic topo- 
logically ordered phases can appear, including f racti onal 
quantum Hall liquids such as the Laughlin stat e 1 6 * 10 ! and 
more exotic Hall liquids which rely on the presence of the 
lattice™] 

A related class of lattice-models relies on spin-orbit 
coupling to generate complex hopping terms in single 
particle tight-binding Hamiltonians. If the resulting 
single-particle bands are flat and have a non-zero Chern 
number, these systems can support states resembling 
fractional quantum Hall (FQH) liquids that are known as 
fractional Chern insulators (FCI) QUEI Fractional Chern 
insulators have been most convincingly shown to exhibit 
the same type of topological order as FQH states by 
analysing their particle entanglement spectraP^ 

The entanglement spectra (ES) have been initially in- 



troduced by Li and Haldane 19 in the context of t he 
FQHE, stimulating an extensive range of studies P^OH-HI 
They have also been studied and applied to several other 
systems includin g sp in systemsp^HSl as well as topolog- 
ical insulatorsf 46 48 Bose-Hubbard models 49 or complex 
paired supcrfluids. 50 The ES corresponds to the spectrum 
of the reduced density matrix of the system groundstate 
when one cuts the system into two parts. The system 
partition can be performed in different manners such as 
a real space, momentum or particle space partition. Each 
cut can unveil different aspects of the state that is probed. 
In the case of the FQHE, the ES are related to the bulk or 
edge excitations. As these features characterize the given 
phase, the ES acts as a fingerprint of the system that only 
requires knowledge of the groundstate wavefunction. 

For the above reasons, the ES were found to be partic- 
ularly well suited as a tool to characterize FCI states as 
they only require knowledge of the ground state wave 
function which can be obtained numerically for small 
model systems.^ Hence, it was possible to establish a 
detailed correspondence between the entanglement spec- 
tra of FCI with those of fractional quantum Hall states.^ 
In particular, it has recently been showir^l that the ES 
is able to distinguish between a Laughlin-like state and 
a charge-density wave state (CDW). However, an im- 
portant difference is that for fractional quantum Hall 
states, very accurate analytic many-body trial wavefunc- 
tions capturing the essence of these strongly correlated 
quantum liquids are knownJS2HH] so a very detailed un- 
derstanding of the topological order and the fundamen- 
tal excitations in FQH systems has been achieved. We 
note that mappings of FQH wavefunctions onto topo- 
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logical flat bandipSiSll have recently led to some encour- 
aging resu lts, in cluding considerable overlaps with FCI 
eigenstate d 59 * 60 ! and the demonstr ation of an analytic 
continuation between these systems! 60 * 62 ! 

In this paper, we focus on quantum Hall states on lat- 
tices with a homogeneous density of gauge flux, as this 
gives us access to a lattice based system where quantum 
Hall states can be understood both in terms of entan- 
glement spectra and many-body trial wavef unctions. As 
opposed to FCIs, quantum Hall states on lattices ad- 
mit a continuum limit which is known to be exactly the 
usual FQHE. We deploy exact numerical diagonalization 
on a square lattice with periodic boundary conditions 
as the main tool for our investigation. Our study pro- 
vides an independent identification of the incompressible 
phases on the lattice as fractional quantum Hall states 
by analysing the counting of quasiparticle states in the 
particle entanglement spectrum. To obtain the expected 
state counting, we establish the correspondence between 
the momentum sectors of a lattice-based system with pe- 
riodic boundary conditions and the corresponding con- 
tinuum problem on the torus. We find that for small 
enough particle density n per lattice site, the counting 
of the continuum problem is accurately reproduced. In 
particular, this allows us to mount evidence in favour 
of the non-Abelian Moore-Read state as the groundstate 
of bosons with repulsive three-body interactions at fill- 
ing factor v = 1 on the lattice. Finally, we study how 
the entanglement gap, defined as the distance between 
the eigenvalues of the entanglement spectrum related to 
the universal property of the bulk excitations and the 
non-universal states at higher entanglement energy, re- 
lates to the magnitude of the overlap with the respective 
trial wavefunctions describing the target phase. This es- 
tablishes the entanglement gap as a good proxy for the 
overlap, and hence for the stability of topological order. 

The structure of this paper is as follows: In section |H} 
we introduce the model Hamiltonian for bosons on lat- 
tices in the presence of magnetic fields and with periodic 



boundary conditions. Section III is devoted to a discus- 
sion of the corresponding continuum FQHE problem on 
the torus, where we introduce the many-body wavefunc- 
tions of the target phases which we explore in this paper, 
namely the Laughlin state (section III A I and the Moore- 



Read state (section III B ) . We then review the definition 
and characteristics of entanglement spectra in section [TV] 
and present a detailed analysis of the particle entangle- 
ment spectra of our two target phases in section |Vj In 
the case of the Laughlin state, we describe how certain 
features of the ES can be used as a probe to detect a 
competing Bose-Einstein condensed phase. Finally, our 
results are summarized in section IVTl 



II. MODEL 

We study the physics of interacting bosons on a two- 
dimensional square lattice in a homogeneous magnetic 



field applied in the direction perpendicular to the lat- 
tice. This problem is described by the Bose-Hubbard 
model with minimal coupling to a gauge field by Peierls' 
substitution. We further assume the presence of on-site 
two-body interactions of strength U and three-body in- 
teractions of strength V, yielding the many-body Hamil- 
tonian 



U = -t J] {e iA **'a\a Tl +h.c.) 

<r,r') 

+ U alala r a r + V alalala T a r a T . 



(1) 



Here, (r, r') denotes neighbouring lattice sites r = (x,y), 

&r are bosonic annihilation (creation) operators, and 

A rr i = A ■ dl are Aharonov-Bohm phases deriving 
from the coupling to the underlying vector potential A. 
We adopt units such that the lattice spacing is one, and 
positions (x,y) can be indicated as integers. 

Experimentally, bosonic Hubbard models can be engi- 
neered in optical lattices systems,^ where gauge poten- 
tials can be simulated by a range of different setups 
Two-body interactions can b e conveniently introduced 
by Feshbach resonances ] 64 * 65 ! and there are proposals for 
three- body interactions based on strong three-particle 
losses! 66 67 Given this pace of progress in simulating Hub- 
bard models, Hamiltonians of the form ([lj may be real- 
izable within the near future. 

For our numerical exact diagonalization calculations, 
we express the Hamiltonian ([!]) for an ensemble of N 
bosons on a square lattice with N s = L x L y sites in the 
presence of N$ flux quanta and with periodic boundary 
conditions in both the x- and y-direction. This set-up 
corresponds to a field of flux density — N^/N s , which 
we choose to describe in the Landau gauge 



A = lix-n^xe. 



(fiXHy , 



(2) 



such that momentum in the y-direction is a conserved 
quantity. Due to the concurrence of a periodic lattice 
potential, periodic boundary conditions, and the pres- 
ence of a magnetic field, this translational symmetry 
is reduced^ to possible momenta of k y = 2rrm/K™ ax , 
m = 0, . . . , _ftT™ ax — 1 with the maximal momentum in- 
dex given by the greatest common denominator 



^max 



gcd(A^,L a 



(3) 



This can be simply explained by applying Blochs theo- 
rem to magnetic unit cells enclosing an integer number 
of flux quanta. Due to the reduced symmetry, orbitals 
are labelled by a sublattice index s for the y-position in- 
side the magnetic unit cell in addition to the momentum 
k y . Both interaction terms in ([T]) conserve this Landau- 
momentum, so the Hamiltonian is block-diagonal and 
we construct the eigenstates in the Fock space given by 
\ a ) = Y[^\) n<; ^W) ^ where C = {x,k y ,s) denotes the 
set of single-particle quantum numbers. 
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We should stress that this set of states does not im- 
ply a projection to the lowest Landau-level. Instead, 
it includes all of the bands of the fractal single-particle 
spectrum known as the Hofstadter butterflypO such that 
Landau-level mixing is part of the model. The equivalent 
Landau level filling, i.e. the particle density with respect 
to the number of states in the lowest band, is given by 
v = N/N <l> . 



III. FQH ON THE TORUS 

We now briefly describe some properties of the FQHE 
in the torus geometry. We consider a torus spanned by 
Li = L' x e x and L 2 = L y e y , where e x and e y are two per- 
pendicular unit vectors. Given a setting where the torus 
is pierced by A^ flux quanta, we have L' x L' y = 2 / jrl B N$, 
where Ib is the magnetic length. The Hamiltonian is 
given by 



+V % - ~ r*) (4) 

i<j<k 

where Hi = —ih\Ji — eA(r,) is the canonical momentum 
of particle j in the presence of a magnetic field. Since 
we use periodic boundary conditions, the delta function 
is defined as <5(r) = J2 n , m ^( r + + TO L 2 ). 

In the Landau gauge (A = 2TTn^xe y ), the one-particle 
orbital in the lowest Landau level with momentum index 
j = 0. . . Na, — 1 is given by 



1 



exp 



(y - ix) 



21% 



where $ 



«7rr(n+a) +227r(n+a) (z-\-b) 



(5) 



are 



the generalized Jacobi thcta functions. This model is 
the continuous version of 0. 

As for the lattice, the momentum in the y-direction 
is a conserved quantity and the A^-body Hamiltonian is 
block diagonal with respect to the total momentum IC y = 
ii modAT^. Note that in the case of the lattice, the 
total momentum JC y is defined modulo AT™ ax . 



A. Laughlin state 

In the lowest Landau level, when only two-body in- 
teractions are present i.e. V = in equation Q, the 
Laughlin state is the densest zero energy groundstateP^ 
A hallmark of this phase is its twofold groundstate de- 
generacy. For a finite size system with A^ = 2N the two 



0, N. 



groundstates can be found at momenta of /C, 
They are given by: 

*(zi, ...,z N ) = / r eiOi, • ■ • , z N )F c . m .{Z)e-^* x t /l z , 

(6) 

where F c m is a center of mass wavefunctions that de- 
pends only on the center of mass coordinate Z = Zi, 
while / ro i is the wavefunction describing the relative mo- 
tion. On a rectangular torus of size (L' x x L' y ), we have 



i<j 



L' 



(7) 



Due to the symmetry under translations of the center 
of mass, the center of mass wavefunction at v = 1/2 is 
two-fold degenerate^ and is given by 



L i N *~ 2 

2 ' 4 
2-iVj, 



2Z 
Li, 



2i- 



■L' x 

>,— 
U, 



(8) 



where / = 0, 1 indexes the two degenerate wavefunctions. 



B. Moore- Read state 

In the lowest Landau level, the Moore-Read state^ 
is the densest zero-energy ground state of the hardcore 
three-body interactions, given by Eq. Q with U = 0. 
It embodies the physics of a chiral p-wave superconduc- 
tor of composite fermionsj 71 * 72 ! which can be cast in the 
real-space form of a BCS paired state in terms of the pair- 
wavefunction l/(zj— Zj). An equivalent expression can be 
found for the torus, where the ground state is three fold- 
degenerate. On a rectangular torus, these groundstates 



can be found at fc y -momenta K. y 



{0,0, f }. For our 



purposes, it is most useful to obtain the Moore- Read trial 
wavefunctions from the Laughlin state by using the Cap- 
pclli formulapS which relates the Moore- Read state of N 
particles to two independent layers of Laughlin v = 1/2 
states with half the number of particles. Using notations 
for a sphere or disk: 

* Pf (zi, . . .,z N ) = 

I N/2 

s n (^-^) 2 (^ +l -^+,) 2 1 • (9) 

\i<j=2 

where S is the symmetrization operator. On the torus, 
the Laughlin state is two fold-degenerate while the 
Moore-Read state is three fold-degenerate. Similar to 
what happens for the quasihole states on the disk or 
sphere geometry, the symmetrization induces linear de- 
pendencies. Here, one can write four states: 



V% T (z u ...,z N ) 



S(^ K 9 T (z 1 ,...,z N/2 ) 

VI 

X*^t Oat/2+1, ■•■ 

V2 



(10) 



,zn)) 
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with Ky ± ,Ky 2 being equal to one of the degenerate 
Laughlin states with = or = N' = f . 

The total momentum Ky = K y v i + K^modN^ per- 
fectly matches the one of the Moore-Read state: taking 
K^ x = K^ = yields a first K£ = MR state, taking 
= K^ 2 = N' again yields K? = and finally tak- 
ing Ky x — and Ky 2 = N' one obtains the remaining 
MR ground state at = N' = N/2 state (Choosing 
Ky t — N' and Ky 2 — is equivalent due to the sym- 
metrization and does not yield an additional state). 



An equivalent classification can be chosen for the contin- 
uum fractional quantum Hall problem on the torusp^ ex- 
cept that the corresponding Landau-momentum ICy can 
take all iV^ distinct values. To compare the entangle- 
ment spectra in these two distinct situations, we identify 
the momentum sectors modulo if™ ax . In particular, we 
will study the counting of the number J\f of low-lying en- 
tanglement eigenvalues below an entanglement gap. In 
this case, we expect the following mapping between the 
values for the torus Nt and lattice Ml 



N L {k y ) 



3d K" 



)N T {K2 



(12) 



IV. PARTICLE ENTANGLEMENT SPECTRA 

The entanglement spectrum gives access to many of 
the spectral properties of the system which are encoded 
in the groundstate wavefunctionP^ It is defined from the 
reduced density matrix of a subsystem resulting from 
the partition of the system into two (or more) parts A 
and B. For the particle entanglement spectrum (PES) 
this partition consists in distributing the particles into 
two subgroups (A and B) while keeping the geometry 
unchanged.^ The reduced density matrix pA = Tr^p, 
obtained by tracing out the Ng particles in the B parti- 
tion, yields the entanglement spectrum by diagonalising 
and classifying the resulting eigenstates according to the 
symmetries of the problem. This process is equivalent 
to a Schmidt decomposition of the original many-body 
state into orthogonal bases for the partitions 



!*> = ££< 



-£w,i/2 



I* 



(11) 



where w stands for quantum numbers designating a sec- 
tor of the decomposition and i indexes states in each 



sector, and the eigenvalues X z 



S ro ,i/2 f de- 



composition ar e re presented on a logarithmic scale. The 
vectors in Eq. jllj ) are orthonormal i.e. iV&^j' j) — 
(vf^Jv!^/ j) — Si.jS-^^'. The entanglement spectrum is 
given by plotting £'s over the relevant w. 

It has been observed that model states such as the 
Laughlin or MR states have a characteristic PES*^ the 
number of non-zero eigenvalues for pa is identical to the 
number of quasihole states for a similar system with the 
same geometry and Na particles. This number is usu- 
ally exponentially lower than the dimension of pa- The 
same features can often persist for eigenstates of realis- 
tic interaction Hamiltonians: for groundstate wavefunc- 
tions with robust topological order, one should observe a 
clearly defined entanglement gap, separating an ensemble 
of low-lying entanglement eigenvalues from non-universal 
eigenvalues located at higher entanglement energies. No- 
tice that model wavefunctions can be thought of as hav- 
ing an infinite entanglement gap. 

For the model associated to the Hamiltonian |IJ), we 
classify the sectors of the entanglement spectrum by the 



Note that for degenerate groundstate manifolds, the 
PES has to be calculated for the incoherent average re- 
duced density matrix p to t for the ensemble of ground- 
states {I^q)}, given by the sum 



Ptot 



des 



(13) 



over the das degenerate ground states. As discussed in 
Ref. [7H this definition yields model state PES for degen- 
erate groundstate manifolds which recover the properties 
of the non-degenerate case on simply connected surfaces. 



V. TARGET PHASES 

The presence of incompressible fractional quantum 
Hall liquids is well established for the Hamiltonian ([Tl. 
These states include the fr actio nal quantum Hall liquids 
of the continuum problem] 6 * 10 ! however, the presence of 
the lattice potential also gives rise to additional incom- 
pressible statesP^E^In this paper, our aim is to establish 
the use of entanglement spectra for lattice based systems, 
so we shall focus on the states with an equivalent in the 
continuum case and undertake a comparison of their fea- 
tures. 



A. Laughlin State 

We begin our analysis with the Laughlin state of 
bosons at v = 1/2, as the best investigated quantum Hall 
state on lattices. We use the analytic form of the Laugh- 
lin states in the continuum (Eq. [6| and substitute the dis- 
crete lattice coordinates, such that z — a/£o(me x + ine y ) 
for lattice site i = (m, n) . By virtue of the folding of mo- 



i.e. they may either remain at dif- 
ferent fc-points or are both mapped to zero momentum 



menta ([3]), the two ground states QSJ) now occur at k y 
and k y = NmodK™ 3 * 



Landau-momentum k v 



0. 



jy'max 



1, i.e., w 



if N mod Ky iax = 0. In our simulations, we find that for 
lattice systems with sufficiently small particle density, the 
momenta of the numerically obtained ground states are 
in agreement with this prediction. In particular, there is 
an extended regime where a two-fold degenerate ground- 
state with a finite gap is found.^ Hence, we can analyze 
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the particle entanglement spectrum of the groundstate 
manifold according to the total density matrix (13). 



We begin to illustrate the mapping of the entangle- 
ment spectrum on the torus to the lattice (12) for a 
small model system with N — 6 particles moving in the 
field of = 12 flux quanta. We first consider the PES 
with Na — 3. Fig. [T^ shows the PES for the Laughlin 
state on the torus geometry The PES yields the following 
counting for the 12 distinct K, y -momentum sectors on the 
torus: (10,9,9,10,9,9,10,9,9,10,9,9). Fig. [T]d displays 
the PES for the groundstate of the model Hamiltonian 
that gives rise to the Laughlin state but with an addi- 
tional small contribution from a longer range interaction. 
In that case, the PES exhibits a entanglement gap. The 
counting below the gap exactly matches the one of the 
Laughlin state. 

We now consider the lattice model. The counting of the 
Laughlin on the torus is reproduced exactly on a lattice 
of geometry L x ~ 4 and L y = 12 for U/t = 1 and V = 
[see Fig. [2j»] as this lattice retains K™ ax = N^ = 12. 
Notice that the total number of states above the gap 
per momentum sector is much higher in the lattice case 
than in the continuum model shown in Fig. [TJa. This is 
a consequence of the Hilbert space dimension being set 
by N and L x x L y in the lattice case, and by N and 
in the continuum case. So, our confirmation of a clear 
entanglement gap for such a small number of particles in 
the lattice model is even more remarkable. 

For lattices with other aspect ratios, a folding of the 
momentum axis often occurs, in which case the maximum 
momentum is reduced. Several examples for such entan- 
glement spectra are shown in Fig. [2^b-d). For example, 
in panel (b), for the geometry with L x = 7, L y = 8, and 
-f^max = 4 the use of equation ( 12 ) predicts the counting 
(28,28,28,28), which is indeed reproduced. The same 
result is also obtained for the aspect ratio of 12 x 4 sites 
in panel (c). However, we do not always obtain a PES 
with a well defined entanglement gap. For the lattice ge- 
ometry of L x — 6 and L y = 8, we find that no threshold 
value £ t for the entanglement energy yields a clear-cut 
definition of N^ t {k y ). We could speculate whether this 
is due to the commensurability of the number of parti- 
cles with L x . In such geometries, it has previously been 
found that CDW states can intervene 1^21 However, there 
is a range of phases which may be competing with frac- 
tional quantum Hall liquids in optical flux lattices, which 
include bosonic condensates with symm etry breaking,^ 
or more general supersolid phases ES1Z6J 

As a first step towards understanding the lattice which 
does not conform to the picture of an incompressible 
Laughlin state (L x — 6, L y — 8), we investigate sev- 
eral additional entanglement spectra for this system, 
analysing the dependency on the number of particles in 
the partition A. The results for Na — 1 and Na = 2 
are shown in Fig. [3] Firstly, we note that the entangle- 
ment spectrum for Na = 1 carries non-trivial information 
for the lattice, while the corresponding continuum limit 
would yield a number of eigenstates which is given by the 



total Hilbert space dimension, i.e., the number of states 
in the lowest Landau level. On the lattice, we instead find 
that a gap opens in the entanglement spectrum above a 
number of low-lying states which precisely matches the 
number of eigenstates in the lowest Landau level. Specif- 
ically, Fig. |3ja) reveals precisely 12 states below the gap 
located at about £ = 6. For a Laughlin state, we would 
expect all twelve of these states to be degenerate. Sec- 
ondly, we find that there are two eigenstates which are 
separated from the other 10 by a further entanglement 
gap located near £ = 2, reproducing the same number 
and momentum sectors K y = and K y — 2 of the 
two low-lying states that we had observed in the PES 
for Na = 3 in Fig. [2jcL) . To complete our survey, we 
also examine the entanglement spectrum for Na — 2 in 
Fig. j3^b) . Again, we find two degenerate low- lying eigen- 
states with an entanglement gap near £ = 2 and located 
in the same momentum sectors. This invariance of the 
number of entanglement eigenvalues with the number of 
particles in the partition is fundamentally different from 
the behaviour that we expect from topologically ordered 
phases. By contrast, the ability to absorb further parti- 
cles without any change of the physical properties (i.e., 
the number of low-lying excitations) can be seen as an 
indication of the physics of Bose condensation. 

To probe for the presence of a Bose condensate, we use 



calculated 



the single particle density matrix • 
between lattice sites i, j. This matrix is exactly a reduced 
density matrix for the specific value Na = 1. A state 
with a finite condensate fraction is signalled by single a 
large eigenvalue Ao of p s , whose magnitude scales with 
the system size N. However, we find that the groundstate 
wavefunction for our system with L x — 6 and L y = 8 and 
N^ — 12 yields a twofold degenerate pair of largest eigen- 
values Ao = Ai ~ 2.472. This characteristic is known to 
be associated to discrete symmetry breaking in the ther- 
modynamic limit We thus follow the procedure intro- 
duced in Ref. [75] and calculate the density matrix for 
symmetry-broken states which are obtained in our case 
by constructing superpositions formed of the two low- 
est lying eigenstates \S) = c |* ) + (1 - |co| 2 )|*i > that 
optimize the largest eigenvalue of . For the symmetry- 
broken state that results from superposing two states 
with different momenta, we find a single large density ma- 
trix eigenvalue Aq(S') = 4.4911, corresponding to a con- 
densed fraction of 74.8% for the N = 6 particle system 
(the condensate fraction rises to 95% as interactions are 
reduced to U — O.lt). At the same time, the state breaks 
translational invariance, forming stripes running around 
the short cycle of the simulation cell. Similar finite size 
effect had previously been reported for lattices™ as well 
as for continuum problems of bosons.^ It is likely that 
properties of the particular lattice size which we discuss 
here are related to its flux density of precisely = 1/4. 
At this value of n^, the single particle Hofstadter spec- 
trum consists of a single, moderately wide band which 
naturally supports Bose condensation at low interaction 
strength. As we do not examine this question in further 
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FIG. 1. Particle entanglement spectra for N = 6 bosons on the torus of unity aspect ratio at filling factor v = N/N$ = 1/2, 
for a particle partition with Na = 3. Left: bosons interacts through hardcore interaction. Right: bosons interacts through 
hardcore interaction and an additional longer range interaction. 




FIG. 2. Particle entanglement spectra for N = 6 bosons on a lattice at filling factor v = N/N^ — 1/2, for a particle partition 
with Na = 3 and for different lattice geometries. The spectra are calculated for the two-fold degenerate groundstate manifold 
of the Hamiltonian with U/t = 1 and V = 0. 



detail, we can only speculate whether all of the above fea- 
tures survive in the thermodynamic limit, in which case 
the phase could be considered a supersolid. 75 For the 
purpose of the current paper, we can conclude that an 
entanglement spectrum with few low-lying states, whose 
number remains invariant for different partitions of the 
system, is indicative of a condensed state. Cases where 



such eigenvalues occur in different momentum sectors are 
likely related to condensates with symmetry breaking. 

Returning to our main discussion of the properties of 
the Laughlin state, we present a collection of the prop- 
erties of the entanglement spectra in table [TJ which gives 
an overview for several lattice geometries that we have 
studied. To summarize our principal findings from these 
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FIG. 3. Particle entanglement spectra for N = 6 bosons on a lattice of 6 x 8 sites at filling factor v = N/N$ = 1/2 for particle 
partitions with Na = 1 (a) and Na = 2 (b) [Na ~ 3 shown in Fig. [2^d)] . For Na = 1, the upper blue dashed line marks the 
entanglement gap below which the counting matches the lowest Landau level state counting of 12 states, and the lower black 
dashed line marks the entanglement gap that isolated the two low-lying eigenstates. For Na = 2, the black dashed line marks 
the entanglement gap that isolated the two low-lying eigenstates. 
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TABLE I. Properties of the particle entanglement spectra of the groundstate manifold of the Hamiltonian for U/t = 1 and 
V — |T| for different model systems of N bosons with N$ = 2N flux quanta. We indicate the momenta of the two 
degenerate groundstates fc^GS) as well as the energy gap A and the groundstate energy splitting 5. The counting of the 
particle entanglement spectrum is shown for the partition with Na = [N/2\ and for different lattice geometries. For both these 
properties, we indicate the agreement with the predictions for the Laughlin state and the entanglement gap Af . 



data, we have established that the counting of excitations 
encoded in the PES for the Laughlin states on the lattice 
Hamiltonian at filling factor v = 1/2 agrees well with 
the data for the continuum problem in the lowest Lan- 
dau level for small flux/particle density per plaquette. In 
particular, the entanglement spectra of this state show a 
clear entanglement gap Ag. 

Unlike the problem on the torus, where it is custom- 



ary to consider the projection of the Hamiltonian into 
the lowest Landau level, the full lattice Hamiltonian 
includes all Landau (Hofstadter) bands. Hence, we can 
study the effect of band mixing that occurs as a function 
of the interaction strength U. The evolution of the en- 
tanglement gap Aj with U is shown in Fig. HI alongside 
the energy gap A. Unlike the energy gap which always 
increases with U, the entanglement gap reaches a maxi- 



N = 5, L. = 6, L, = 1 
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0.15 
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entanglement gap 




FIG. 4. Energy gap A and entanglement gap A^ in the Na = 
2 and in the Na = 3 sectors as a function of U for v = 1/2, 
AT = 5, L x — 6 and L y = 10. The vertical purple line is the 
band gap. Note the different offsets on the scales for A and 



FIG. 6. Entanglement gap A^ as a function of flux density 
for the Laughlin state with N = 5 particles. The vertical 
dotted line indicates the critical value n£ ~ 0.4 up to which 
the groundstate exhibits the Chern number of the Laughlin 
stateP 



mum value for an interaction strength of the order of the 
band gap and then decreases. 
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FIG. 5. Overlap O and entanglement gap A^ in the Na = 3 
sectors as a function of U for v = 1/2, N = 5, L x = 6 and 
L y = 10. The vertical purple line is the band gap. Note the 
different offsets on the scales for O and A^. 

We also computed the total overlaps to t = 
3 Si,j=i I (*GS,i|*modei,j) | 2 of the exact groundstates 
with the model state as a function of the interaction 
strength. The results are shown on Fig. [5j One can 
notice that the overlaps are very high. Moreover, the 
overlap is an increasing function of U, as is the energy 
gap. 

It is now well established that the presence of an en- 
tanglement gap, in conjunction with the specific state 
counting in the PES, characterizes the topological order 
in the system. We are therefore interested to test how 
this measure compares to other signatures of topological 
order, such as the presence of non-zero Chern number for 
the groundstate manifold. A prior study of the Laughlin 
state on lattices had shown that the combined twofold 
ground state manifold has a Chern number of 1, or 1/2 



per state, up to a critical flux density of ~ 0.4.^ We 
now study how the entanglement gap varies as the flux 
density changes, by calculating entanglement spectra for 
systems of constant N on lattices of different geometries. 
The results, shown in Fig. [6j show a full agreement with 
Ref. EH For large U, we find that the entanglement gap 
also closes at ~ 0.4. However, while the Chern num- 
ber jumps instantaneously between integers, the entan- 
glement gap can capture how the topological protection 
of the Laughlin state is gradually weakened and finally 
collapses. 



B. Moore-Read State 

We next considered the Moore- Read (MR) state at v = 
1. In the continuum, as explained in |III B[ the three 
MR states can be obtained from the Laughlin state by 
symmctrization. On the lattice, the same scheme applies 
except that momentums are now defined modulo X™ ax , 
and we use the Laughlin states on the lattice as defined 
in Eqns. ([6|-([8| as the starting point. 

Our numerical work on the lattice is based on the exact 
diagonalization of the Hamiltonian ([I]) , using three-body 
contact interaction which are analogous to the contin- 
uum case, i.e. choosing U = and V/t = 1 in equation 
([!]). Given these parameters, we generally found that the 
groundstate is approximately three-fold degenerate and 
the sectors in which the three ground states appear are 
given by the expected momenta, subject to the folding 
rule (12 1. For geometries where the latter is satisfied, 



we compute the particle entanglement spectrum of the 
groundstates total density matrix. As in the Laughlin 
case, the particle entanglement spectrum is gapped and 
the number of states below the gap is given by the one 
predicted from the folding rule and the torus counting 
[see Fig. [7ja)]. The results for the different systems we 
studied are gathered in table [Til 
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TABLE II. Properties of the particle entanglement spectra of the groundstate manifold of the Hamiltonian for (7 = and 
V/t = lffl for different model systems of N bosons with N$ = N flux quanta. We indicate the momenta of the three 
degenerate groundstates k y (GS) as well as the energy gap A and the groundstate energy splitting S. The counting of the 
particle entanglement spectrum is shown for the partition with Na = [N/2\ and for different lattice geometries. For both these 
properties, we indicate the agreement with the predictions for the Moore- Read state and the entanglement gap Aj. 



These results show that it is theoretically possible to 
obtain the Moore-Read phase on the lattice using three- 
body contact interactions. Even though three-body in- 
teractions can be realized for cold atoms using more elab- 
orate experimental settings^ 66 * 67 ! the most relevant inter- 
action for bosons on a lattice is the two-body hardcore 
interaction. Thus, we wonder if the Moore-Read state 
can also be stabilized with this type of interaction. In 
the continuum limit and in the lowest Landau level ap- 
proximation, they are several numerical eveidences that 
such a phase can be stabilizePsHUl Given the presence 
of LL mixing mixing in our model, we also establish 
to which extent this mixing affects the stability of the 
Moore-Read state. To answer these questions, we diag- 
onalize the Hamiltonian ([I]) with V = at v = 1. For 
small interaction strength U, the energy spectrum ex- 
hibits the correct groundstate degeneracy and -sectors. 
In addition, the entanglement spectra exhibit the correct 
state counting, albeit with a smaller entanglement gap 
than in the case of three-body interactions. A closer look 
at the spectrum reveals, as displayed in Fig. [HJ that the 
energy gap and the groundstate energy splitting are of 
the same order: While a fully realized MR phase should 
have a very small spread to gap ratio, this is clearly not 
the case here. 



We find that the gap closes for U c ~ 1.25t, while the 
spread between the ground states grows monotonically 
with U. Hence, the spread to gap ratio is a rapidly in- 
creasing function of U indicating how sensitive the topo- 
logical degeneracy is to the strength of two-body inter- 
actions. At the flux density = 1/6 shown in Fig. |8j 
the lowest Landau level is still very narrowp^H so we 



N = 6, L x = 8, Ly = 6, N A = 3 




-1 1 2 3 4 5 6 

FIG. 7. Particle entanglement spectra for N = 6 bosons on a 
lattice at filling factor v = N/Nj, = 1, for a particle partition 
with Na = 3 and for different lattice geometries. The spectra 
are calculated for the three-fold degenerate groundstate man- 
ifold of the Hamiltonian ^ with [7 = and V/t = 1. Left: 
The counting of the states below the gap is (7,6,6,7,6,6). 
It matches the one of MR quasiholes states on the torus as 
K™* = N*. 



interpret this strong suppression of the gap as result- 
ing from Landau level mixing. Furthermore, the single 
particle gap between the LLL and second LL is about 
A sp ~ 1.68i at n.0 = 1/6, which is of the same order as 
the energy scale of interactions nU c at the point U c where 
the gap closes. We can conclude, at least for small sys- 
tems, that the Landau level mixing resulting from large 
two-body contact interactions tends to destroy the MR 
phase. 

Finally, we have computed the entanglement gap of the 
total matrix density, taking the two lowest energy state 
in the K y — sector and the lowest in the K y = N/2 
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FIG. 8. Energy gap A = E3 — E2 and groundstate energy 
splitting 5 = E2 — Eo as a function of U for v = 1, N = 6, 



6 and L v 



as well as their dimensionless ratio. 




tanglement spectrum. These systems provide a well con- 
trolled environment away from the pure model states, 
which allows a better understanding of the properties 
of particle entanglement spectra (PES) in quantum Hall 
systems. We have focused on the filling factors v = I and 
v = 1 where the Laughlin state and the Moore- Read state 
should respectively emerge. In both cases, the PES was 
able to discriminate the nature of the state. This result is 
even more remarkable given that the size of the Hilbert 
space, set by the number of particles and lattice sites 
(rather than flux quanta), is exponentially larger than in 
the continuum limit. Interestingly, the PES was able to 
give insight about a competing Bose-Einstein condensate 
phase, which we have associated with low lying entangle- 
ment eigenstates whose number is invariant under the 
number of particles in the partition. We have also shown 
that the entanglement gap collapse in the PES predicts 
a critical density of flux below which the Laughlin's 
physics emerges; our value of is in agreement with a 
previous study based on Chern numbers. We have used 
the PES to confirm the realization of a Moore- Read state 
at v = 1 in the presence of on-site two-body contact in- 
teractions, only. Furthermore, we have given evidence of 
how Landau level mixing arising from these two-body 
contact interactions tends to destroy the bosonic MR 
state as its magnitude is increased. 



FIG. 9. Overlap O and entanglement gap in the Na = 3 
sectors as a function of U for v = 1, N = 6, L x = 6 and 
L y — 6. Note the different offsets on the axes for the overlap 
O and entanglement gap A^. 



sector even when they were not the three lowest energy 
eigenstates. For reference, we also evaluate the overlap 
of these states with the model MR states, as discussed 
above. The results, shown in Fig. [9j are consistent with 
the previous conclusions: The phase obtained at small 
U is most likely the MR phase, with large overlaps at 
small U, but the phase is destroyed by increasing the 
interaction strength. 



VI. CONCLUSION 

In this paper we have analyzed the bosonic fractional 
quantum Hall states on lattices through the particle en- 
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